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ABSTRACT 

We give an example in which it is possible to understand quantum statistics using 
classical concepts. This is done by studying the interaction of charged matter oscilla- 
tors with the thermal and zeropoint electromagnetic fields characteristic of quantum 
electrodynamics and classical stochastic electrodynamics. Planck's formula for the 
spectral distribution and the elements of energy huj are interpreted without resorting 
to discontinuities. We also show the aspects in which our model calculation complement 
other derivations of blackbody radiation spectrum without quantum assumptions. 
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Recently, Tersoff and Bayert-'-] have shown that the statistics obeyed by the particles 
cannot be used as a criterion for their distinguishability. Here we show that the general 
proposal of Tersoff and Bayer, when adapted to the blackbody radiation phenomenon, 
provides a clue for bringing quantum and classical physics into a closer relation. In 
this regard, the similarity between quantum electrodynamics (QED) and the classical 
theory called stochastic electrodynamics (SED) has been stressed by several authors in 
recent years P^^] . It should be mentioned that to obtain a closer relation between the 
classical and quantum approaches was one of the goals of Planck and Einstein research 
in the first and second decades of this century. It is well known that Einstein himself 
never gave up to obtain a deeper understanding of the light-quantum concept within 
the realm of classical electrodynamics [^1 . Following then the proposal of Tersoff and 
Bayer, we show that it is possible to reconcile the statistics of Bose and Einstein with 
classical electrodynamics. As a natural byproduct of our analysis, we interpret Planck's 
radiation law for the blackbody radiation, without using concepts strange to classical 
physics, like discontinuities in the energy of the charged particle during the emission 
or absorption processes. This reinforces the attempts to interpret the "photon" model 
using only the undulatory aspect of the electromagnetic field, thus providing a simpler 
picture of the "light-quantum" [^] . 

For a single one-dimensional nonrelativistic matter oscillator (mass m , charge 
e and frequency co such that e^cu/mc^ <S 1) interacting with a cavity radiation 
with spectral distribution p{u!,T) — p^iuj) + Pj.{oj) , there is a simple relation (valid 
within both QED and SED) between the average oscillator energy (e) and p{uj,T) , 
namely [■'■^l 
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where a;inax is the maximum frequency compatible with nonrelativistic motion. 

In (1) we have considered that huj/2 is the average zeropoint (or zero temperature) 
energy of the oscillator, (u) is its average thermal energy and 

is the spectral distribution of the zeropoint electromagnetic field (which exists in both 
SED and QED). We are denoting by pj,{oj) the spectral distribution of thermal 
radiation with temperature T . As far as we know, formula (1) was firstly obtained 
by Planck (see ref. [8], p. 870) with p^iuj) = . The inclusion of the zeropoint field 
contribution to the oscillator average energy was made later, in 1911, also by Planckf . 

According to SED, the spectral distribution p{uj,T) which appears in (1), is 
associated to the random fluctuations of the classical electromagnetic fields existing 
inside the cavity. These fields are correlated, and one of the correlation functions is 
defined through the ensemble average 

[3,4] 

{E{t) ■ E{0)) 
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where E{t) is the electric field at time t , in the origin of coordinate system (oscillator 
position). It is also assumed that {E{t)) = 0. We shall adopt the classical stochastic 
electrodynamics approach in what follows. 

Within the realm of this enlarged version of Maxwell's electromagnetic theory it 
is possible to derive the exact expression for p(a;,T) without quantum assumptions. 
This is well known and a number of researchers have published, and commented, several 
(different) derivations many times in the pastC^L To our knowledge the first one was 
published by Einstein and Sternlll] in the paper entitled "Some Arguments for the 
Assumptions of a Molecular Agitation at Absolute Zero" pubhshed in 1913. According 



to Einstein and Stern, "the assumption of the zeropoint energy opens a way for deriving 
Planck's radiation formula without resorting to any discontinuities whatsoever" . 

For completeness we shall outline here another simple classical derivation which 
is based on the interaction of a single oscillator with cavity radiation. According to 
this derivation, which was discussed previously by some authorsC^' ''"^], the starting 
point is the Einstein formula for the thermal fluctuations of the oscillator energy. 
In agreement with the Einstein formula, which may be interpreted as definition of 
temperaturet^"^] (see also ref. [8] p. 874), the total oscillator energy has thermal and 
zeropoint fluctuations. Therefore, it is possible to show that the energy variance is 
such that [13] 

{e')-{ef = {ef = ^T^ + i^)' ' (4) 

where (e) is given by (1). As we already said, this result was used previously by 
several authors (see for instance Boyert-'-'^] or Pranga and Santost^^]). In equation (4), 
the first term of the last equality represents the thermal fiuctuations. The second term 
is associated to the zeropoint fluctuations of the oscillator energy. It is assumed that 
the thermal and zeropoint fluctuations are statistically independent (see de la PenaW , 
p. 475) in the sense that their variances simply a , as is shown in equation (4). 

The differential equation (4) has an exact solution for the average energy (e) of 
the oscillator as a function of the temperature T . This solution is precisely t-'-'^] 

, , huj , f hu \ huj huj . , 

as it is easy to verify. The hmit kT ^ hu gives the expected result, namely 
(e) ~ {u) ~ kT , and the low temperature limit is also correct. Besides these 
results, one can obtain the spectral distribution p{u!, T), by combining (5) and (1). To 



our knowledge this is one of the simplest classical derivations of the Planck spectral 
distribution. An important remark is that the result (5) follows directly from (1) and 
(4), and the equations (1) and (4) follows from the classical stochastic equations of 
motion for the charged oscillator [^~^]. 

In order to make a closer comparison with the original 1900 Planck's paper, and to 
make the connection with the statistics of Planck, Boset^^] and Einstein, it is necessary 
to study the same physical system considered by Planck or a similar one. We may 
consider, for instance, a macroscopic group of A identical matter oscillators, which are 
coupled to each other hke the ions in a sohd lattice. We assume, for instance, that these 
oscillators are in the walls of a cavity which exist in a solid material. The cavity has 
thermal and zeropoint radiation, which are in equilibrium with the matter oscillators 
of the surrounding walls. Therefore, this physical model will be naturally related to 
the Einstein and Debye models for the thermal vibrations of a solid body. 

A classical interpretation of the Debye law for the specific heat of sohds was pre- 
sented by Blanco, Franga and Santos[16] quite recently. These authors have used the 
SED approach, thus presenting a derivation of the Debye specific-heat law, using clas- 
sical equations for Brownian ions coupled by linear forces. Their work "provides a 
classical interpretation of the specific heat as being due to a continuous distribution of 
energies of the normal-modes vibrations" . The spectral distribution of these vibrations 
is given by p{uj, T) = {hw^ coth.(huj/2kT). Therefore, this wor 

k[16] 

shows "that 

the close connection (discovered by Einstein in 1907) between the Planck blackbody 
spectrum and the specific-heat law is maintained even if both are interpreted along 
classical ideas" . We recommend the paper by Blanco et al. 

[16] 

to the interested reader. 

Here, however, we shall consider a simpler system (the Planck original system), 
consisting of A identical uncoupled matter oscillators. That is, we shall assume that the 
oscillators are separated from each other. Consequently, each oscillator interacts only 



with the electromagnetic waves associated to the thermal and the zeropoint radiations. 
In this sense, our model is physically equivalent to a rarefied gas of polarizable molecules 
in equilibrium with blackbody radiation. Analyzing this simple system (uncoupled 
oscillators plus radiation) along the ideas proposed by Tersoff and Bayer, we show 
that it is possible to understand better the relation between Maxwell electromagnetic 
theory, Planck's radiation law and Bose-Einstein statistics. In this sense our paper 
complements the other derivations (classical and quantum) of the blackbody radiation 
spectrum. 

Let us denote the thermal energy of the A oscillators by U. Therefore, it is possible 
to write U as 

= E«i , (6) 

i=l 

where Ui is only the thermal part of the total energy of the i-th oscillator. The index 
i denotes the site, or position, of the oscillator in space. The zeropoint energy of the 
oscillators is not included in (6). 

On physical grounds, in order for the temperature T of a system of A oscillators 
to be meaningful, we expect the system to be in equilibrium, the average total thermal 
energy (U) of the oscillators to be constant and we also expect the system to be very 
large (A 3> 1 or A — > oo). According to (1) and (6) the average thermal energy is 
simply: 

(U) = A{u) , (7) 

with (U) ~ 0, if T ^ 0, or (U) AkT if kT hu . However, the actual energy 
U is not constant. It has fluctuations, due to energy exchange with the zeropoint and 
the thermal radiations, as well as due to the flow of energy in and out of the box, or 
volume, containing the oscillators and the electromagnetic radiation. Therefore, we 



expect U to be a random variable which varies continuously within some interval, 
namely 

^min < ^ < ^max , (8) 

where C/min is close to C/max for a large system (A —>■ oo). Nevertheless, we shall 
assume that t/max > f^min and that t/min —>■ for very low temperatures {T ^ 0). 
Moreover, we expect the fluctuations in the thermal energy U to be more important 
for small T. However, the total thermal and zeropoint energy of the system (oscillators 
plus radiation) is assumed to be conserved, via energy conservation laws discussed in 
the recent article by 

Cole[17] , for instance. 
We need to introduce one important step in order to apply the Tersoff and Bayer's 
method and to connect it with the quantum picture and the early Planck's paper. So, 
similarly to Planck, let us formally divide the total thermal energy U of the matter 
oscillators into fractions q , 

' - ^ . (8) 

N being an arbitrary integer that is (at least) of the same order of magnitude as A. 
This subdivision of the thermal energy U is necessary because we want to count the 
number of different ways one can distribute a continuous amount of thermal energy U 
amongst A oscillators. 

Since the oscillator thermal energies can take arbitrary values rt, such that (6) is 
valid, we can introduce a set of numbers {ui} , not necessarily integers (this is an 
striking departure from the Planck's 1900 paper), such that 

Ui = UiQ , (10) 

and 

N ■ (11) 

i=l 



Of course the possibility of integer numbers {n^} is also compatible with (10) and 
(11). In any case, each rii represents the amount of fractions q which an oscillator 
has. We shall assume that these fractions of energy, q = U/N, are distinguishable 
within the context of Tersoff and Bayer classical approach. Moreover, each fraction 
q varies continuously and has fluctuations (see (8) and (9)) and, therefore, cannot be 
strictly identified with the hypothetical Planck's quantum. 

The approach of Tersoff and Bayer requires the assignment of different (random) 
probability weights cij (for the absorption of thermal energy) to identical oscillators 
(i = 1,2, .. .A). So, in order to use their approach, an essential step is to justify this 
assumption. This can be achieved if we recognize that the probability weights must 
be related to the probability of having some thermal energy available to each oscillator 
in its position in space. Regarding this point, our approach will be similar to that 
proposed by Ehrenfestt"'"^] in 1911. Ehrenfest generalized Boltzmann's concept of "a 
priori probability" and used probability weights which vary with the energy interval. 
In the following we shall present a simple model which makes an explicit use of this 
idea. 

Let us assume that the probability weight is proportional to the energy density 
available in site i, which is the i-th oscillator position. As was explained before (see 
(3)), the total electric field, associated to the thermal and zeropoint radiation, in that 
point is a random vector such that 



Therefore, it is quite natural to assume that the probability weights ai may vary 
according to the following simple model: 



im) = 



and 
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where p{Ei) is the probabihty distribution of the classical, Gaussian, random field, 
characteristic of SED (see Boyert"^] and Milonnit'^] , chapters 2 and 8). In other words, 
according to our model, the probability weights ttj may increase (or decrease) with the 
electromagnetic energy density available in the oscillator position. We shall, therefore, 
assume that each ctj may vary randomly from zero (no energy available) to one, 
namely 

< ai < 1 , (14) 

as is expected on physical grounds. 

At this point it is also important to recall that, according to Tersoff and Bayer, 
"it is arguably more natural (i.e., entails a weaker assumption) to take the probability 
weights (ccj) as arbitrary and random than equal" . We shall also assume that 

E«i = 1 (15) 

i=l 

within our model. This is required since the entire system (oscillators plus radiation) 
is in equilibrium. Intuitively speaking, the flow of electromagnetic energy from one 
oscillator site to the other oscillators sites, is compatible with our model equations 
(13), (14) and (15). Moreover, the explicit knowledge of the functional form of the 
probability weights {ctj} is not necessary in order to apply the Tersoff and Bayer 
method as we shall see in a moment. 

The above considerations will allow us to calculate the probability P{ni} , as- 
sociated to the numbers {rij} of distinguishable fractions (g = U/N), which are 
distributed amongst A distinguishable oscillators according to (9), (10) and (11). 
This probability, as proposed by Tersoff and Bayer, is 



where the ensemble average ( ) is related to the continuous and random variation of 
the probability weights {ccj} . In (16) rij! is the gamma function r{ni + 1) if rii is 
not an integer. 

The exact result for this expression, valid even in the general case in which the 
numbers rii are not integers, is simplyt-'-] 



P{ni} 



N\ {A-l)\ 

{N + A-iy. 



(17) 



independently of the numbers {ui} . This remarkable result of Tersoff and Bayer is 
clearly vahd within the realm of SED. It allows the derivation of Planck's radiation 
law without resorting to any discontinuities in the emission or absorption of radiation 
by the matter oscillators, and to interpret Planck's counting (or statistical) methodC^] 
using only the classical properties of electromagnetic fields and mechanical oscillators. 

As A and N are assumed to be very large numbers {A ^ 1 and ^ 1), 
using the Stirling approximation and the Boltzmann relation between the entropy and 
probability, the entropy per oscillator is given by 



S* = - In P{ni} ^ k 



N\ / N 



N , N 
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(18) 



Using (9), the entropy S* can be expressed in terms of the variable U/qA, namely. 



. U \ , L U\ U U 
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(19) 



We shall call this function S* the probabilistic entropy per oscillator. It will be 
identified with the caloric entropy S in the thermodynamic limit A ^ oo. It should 
be remarked that S has to be a function oicu/T, according to the classical Wien lawi'^]. 

Since S* has to be finite, it is expected that the ratio U/qA will be finite in the 
thermodynamic limit. We shall show that the ratio U/qA will be related to the 



properties of a single oscillator, with thermal energy (e) — — , in accordance with the 
SED calculation (see equation (5)). 

The relation between the temperature T and the caloric entropy is well known 
and can be written as 1/T — dS/d{e) . Therefore, according to SED we have 
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where the last equality in (20) was obtained from our previous equation (5). Therefore, 
the simple integration of (20) will lead to 



S ^ k 



(21) 



where (u) = (e) — huj/2. 

Consequently, in the thermodynamic limit A — > oo, one must have 



S* 



s 



(22) 



According to (6), (7) and (8) the ratio qA/U can be written as 



A 



(23) 



because U = A{u) + 5U with {5U) = and \5U\ ~ ^/A{u) for large A. 

Therefore, comparing the expressions (19) and (21) and taking into account (23), 
the limiting process indicated in (22) is equivalent to 



or 
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In other words, hcu is the average value, {q), of the energy fractions used here in 
order to apply the Tersoff and Bayer counting method. As we have already mentioned, 
these fractions q cannot be strictly identified with the hypothetical Planck's quantum 
huj. The result (25) for g is a consequence of the counting method of Tersoff and 
Bayer applied to the thermal energy U which varies continuously and has random 
fluctuations. Moreover, one can show that the value huj obtained for (q) can be 
traced back to the assumed existence of zeropoint electromagnetic fluctuations with 
average energy huj/2 per mode (see (1), (2), (5), (20) and (21)). 

Similar conclusion was obtained by Barranco and Frangat^] in a recent paper, al- 
though through a different reasoning. These authors replaced the Einstein concept of 
"random spontaneous emission" (see also Milonni [4] p. 63) by the concept of "stim- 
ulated emission by the random zeropoint fields" . Using this new concept and the old 
result by Einstein and Ehrenfest t^] , namely, the energy of the molecule (which inter- 
acts with the radiation fields) can vary continuously, Barranco and Pranga extended the 
Einstein (1917) and Einstein- Ehrenfest (1923) works to the realm of classical stochastic 
electrodynamics. As a result, Compton's kinematic relations were interpreted within 
the realm of a classical theory, that is, treating the molecules as classical particles and 
the radiation as classical waves. 

We have seen that, according to SED, the thermal "photons" can be described using 
only the classical aspects of the radiation field. This simplified view of thermal radia- 
tion is important because it is also well known that the photoelectric and the Compton 
effects can be explained without the particle model of the "photon" 

[19] 

. These facts are 

known since the pioneering works by Richardson (1914), Wentzel (1927), Schrodinger 



us 



(1927) and Klein and Nishina (1929) among many othersP^L This suggests to 
that the particle model of the "photon" is not firmly estabhshed and that QED essen- 
tially provides the rules of the interaction without explicitly invoking the corpuscular 
character of the "photon" [^^' (see also Milonnit"^] , pg. 19). It is also important 
to comment another "evidence" of the "corpuscular aspect" of the "photon" (i.e., its 
"indivisibility"). This was offered by the "almost" 100% anticorrelation affecting the 
detection of the "single-photon" (i.e. an experiment with a very weak light beam) by 
two detectors placed on the two output modes of a beam-splittert^-'-]. However, the 
results of this experiment are not conclusive because they may be compatible with 
classical wave models of light [22]. 

Nevertheless, we want to remark that the Tersoff and Bayer derivation of quantum 
statistics based on non-traditional assumptions does not disprove the validity of the 
conventional ones, and that those writers wisely caution that: "Since we have shown 
that either set of counting assumptions can lead to quantum statistics, the statistics 
which "particles" obey cannot be a criterion for their distinguishability" . In this re- 
gard, there are some indications that the distinguishability of "photons" (or classical 
electromagnetic pulses) interacting with optical instruments has been experimentally 
verified. This important discovery, which was reported by at least two different groups 
very recentlyPS] ^ j^g^g j^q^ received the attention it deserves. 

Another remark is that blackbody derivations are notorious as being feasible on a 
number of different models. Therefore, we would advise that a multiplicity of deriva- 
tions from different points of view merely indicates that a successful derivation (con- 
ventional or non-traditional) cannot constitute a proof of any assumed model. In 
this regard we would like to mention other attempts to obtain a classical derivation 
of Planck's radiation law. As far as we know, Einstein and Stern[ll] in 1913, and 



Nernst[24] 

in 1916 were the first physicists who claimed to have derived Planck's for- 
mula without assuming discontinuities (see also ref. [7]). More recently Boyer[^^~^'^] 
and otherst^^^"^"^] have proposed different ways to derive the blackbody radiation spec- 
trum without quantum assumptions. These approaches were based on Planck's second 
theory (1911) which introduced the concept of zeropoint energy of the electromagnetic 
field. It is worthwhile to recall that nowadays these zeropoint fluctuations have been 
extensively observed experimentallyC'^' 34-36] 

According to A. Hermann t"^^] , "Planck realized the radical consequences of his 
formula only much later, and for many years he made attempts to harmonize the 
quantum of action with classical theory. Planck did not believe that he made a complete 
break with the past. He gave the impression that it was impossible to speak of a counted 
probability without a subdivision of energy" . This is in agreement with Ehrenfest old 
opinion that the elements of energy in Planck's calculation were only a formal device: 
"The permutation of the energy symbols {hu) had no more physical significance than 
the permutation of the divider symbols" (see ref. [18] pp. 255 and 256). Ehrenfest's 
formal proof of Planck's combinatorial formula (our equation (17)) was enclosed in a 
letter to Lorentz in 1914. In another (historical) letter addressed to R.W. Wood in 1931 
(see ref. [37] p. 23 for the full reproduction), Planck describes in detail the psychological 
aspects of his approach. "Briefly summarized, what I did can be described as simply 
an act of desperation. I was ready to sacrifice every one of my previous convictions 
about physical laws. Boltzmann had explained how thermodynamic equilibrium is 
established by means of a statistical equilibrium, and if such an approach is applied 
to the equilibrium between matter and radiation, one finds that the continuous loss 
of energy into radiation can be prevented by assuming that energy is forced, at the 
outset, to remain together in certain quanta. This was a purely formal assumption 



and I really did not give it much thought except that, no matter what the cost, I must 
bring about a positive result" . 

According to SED and QED the inevitable loss of energy into radiation is natu- 
rally prevented by the vacuum zeropoint electromagnetic field, which acts as an energy 
reservoir t^"*^] . On the other hand, the Tersoff and Bayer approach, adapted here to 
classical stochastic electrodynamics and the old blackbody radiation problem, also clar- 
ifies (on classical grounds) the statistical method used by Planck. Therefore, Planck's 
formulas (see our equations (5), (17), (19), (20), (21) and (25)) do not necessarily rep- 
resent a complete break with classical physics and the continuity principle. Moreover, 
in their interesting paper, Tersofi^ and Bayer suggest that "quantum mechanics may 
be after all compatible with classical ideas of locahty and distinguishabihty; the key 
appears to lie in a reevaluation of the underlying assumptions about probabilities". 
This is precisely the approach used in SED. In our opinion this general idea deserves 
further attention. 
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